Genetic Analysis of 27 Y-STR Haplotypes in 11 Iranian Ethnic Groups.

BACKGROUND: The study of Y-chromosomal variations provides valuable insights into male susceptibility in certain diseases like cardiovascular disease (CVD). In this study, we analyzed paternal lineage in different Iranian ethnic groups, not only to identify developing medical etiology, but also to pave the way for gender-specific targeted strategies and personalized medicine in medical genetic research studies. METHODS: The diversity of eleven Iranian ethnic groups was studied using 27 Y-chromosomal short tandem repeat (Y-STR) haplotypes from Y-filer® Plus kit. Analysis of molecular variance (AMOVA) based on pair-wise RST along with multidimensional scaling (MDS) calculation and Network phylogenic analysis was employed to quantify the differences between 503 unrelated individuals from each ethnicity. RESULTS: Results from AMOVA calculation confirmed that Gilaks and Azeris showed the largest genetic distance (RST=0.35434); however, Sistanis and Lurs had the smallest considerable genetic distance (RST=0.00483) compared to other ethnicities. Although Azeris had a considerable distance from other ethnicities, they were still close to Turkmens. MDS analysis of ethnic groups gave the indication of lack of similarity between different ethnicities. Besides, network phylogenic analysis demonstrated insignificant clustering between samples. CONCLUSION: The AMOVA analysis results explain that the close distance of Azeris and Turkmens may be the effect of male-dominant expansions across Central Asia that contributed to historical and demographics of populations in the region. Insignificant differences in network analysis could be the consequence of high mutation events that happened in the Y-STR regions over the years. Considering the ethnic group affiliations in medical research, our results provided an understanding and characterization of Iranian male population for future medical and population genetics studies.

SNUPN deficiency causes a recessive muscular dystrophy due to RNA mis-splicing and ECM dysregulation.

SNURPORTIN-1, encoded by SNUPN, plays a central role in the nuclear import of spliceosomal small nuclear ribonucleoproteins. However, its physiological function remains unexplored. In this study, we investigate 18 children from 15 unrelated families who present with atypical muscular dystrophy and neurological defects. Nine hypomorphic SNUPN biallelic variants, predominantly clustered in the last coding exon, are ascertained to segregate with the disease. We demonstrate that mutant SPN1 failed to oligomerize leading to cytoplasmic aggregation in patients' primary fibroblasts and CRISPR/Cas9-mediated mutant cell lines. Additionally, mutant nuclei exhibit defective spliceosomal maturation and breakdown of Cajal bodies. Transcriptome analyses reveal splicing and mRNA expression dysregulation, particularly in sarcolemmal components, causing disruption of cytoskeletal organization in mutant cells and patient muscle tissues. Our findings establish SNUPN deficiency as the genetic etiology of a previously unrecognized subtype of muscular dystrophy and provide robust evidence of the role of SPN1 for muscle homeostasis.

Characterizing Homozygous Variants in Bardet-Biedl Syndrome-Associated Genes Within Iranian Families: Unveiling a Founder Variant in BBS2, c.471G>A.

Bardet-Biedl syndrome (BBS) is a rare inherited ciliopathy disorder characterized by a broad spectrum of clinical symptoms such as retinal dystrophy, obesity, polydactyly, genitourinary and kidney anomalies, learning disability, and hypogonadism. The understanding of the variants involved in BBS-causing genes remains incomplete, highlighting the need for further research to develop a molecular diagnostic strategy for this syndrome. Singleton whole-exome sequencing (WES) was performed on sixteen patients. Our study revealed (1) nine patients carried eight homozygous pathogenic variants with four of them being novel (2) Specifically, a synonymous splicing variant (c.471G > A) in BBS2 gene in six patients with Baloch ethnicity. The identification of runs of homozygosity (ROH) calling was performed using the BCFtools/RoH software on WES data of patients harboring c.471G > A variant. The presence of shared homozygous regions containing the identified variant was confirmed in these patients. In-silico analysis predicted the effect of the c.471G > A variants on BBS2 mRNA splicing. This variant results in disrupted wild-type donor site and intron retention in the mature mRNA. (3) And a deletion of exons 14 to 17 in the BBS1 gene was identified in one patient by Copy-Number Variation (CNV) analysis using the ExomeDepth pipeline. Our results identified the founder variant c.471G > A in the BBS2 gene in the Baloch ethnicity of the Iranian population. This finding can guide the diagnostic approach of this syndrome in future studies.

Clinical application of next generation sequencing for Mendelian disease diagnosis in the Iranian population.

Next-generation sequencing (NGS) has been proven to be one of the most powerful diagnostic tools for rare Mendelian disorders. Several studies on the clinical application of NGS in unselected cohorts of Middle Eastern patients have reported a high diagnostic yield of up to 48%, correlated with a high level of consanguinity in these populations. We evaluated the diagnostic utility of NGS-based testing across different clinical indications in 1436 patients from Iran, representing the first study of its kind in this highly consanguineous population. A total of 1075 exome sequencing and 361 targeted gene panel sequencing were performed over 8 years at a single clinical genetics laboratory, with the majority of cases tested as proband-only (91.6%). The overall diagnostic rate was 46.7%, ranging from 24% in patients with an abnormality of prenatal development to over 67% in patients with an abnormality of the skin. We identified 660 pathogenic or likely pathogenic variants, including 241 novel variants, associated with over 342 known genetic conditions. The highly consanguineous nature of this cohort led to the diagnosis of autosomal recessive disorders in the majority of patients (79.1%) and allowed us to determine the shared carrier status of couples for suspected recessive phenotypes in their deceased child(ren) when direct testing was not possible. We also highlight the observations of recessive inheritance of genes previously associated only with dominant disorders and provide an expanded genotype-phenotype spectrum for multiple less-characterized genes. We present the largest mutational spectrum of known Mendelian disease, including possible founder variants, throughout the Iranian population, which can serve as a unique resource for clinical genomic studies locally and beyond.

Implementation of an In-House Platform for Rapid Screening of SARS-CoV-2 Genome Variations.

BACKGROUND: Global real-time monitoring of SARS-CoV-2 variants is crucial to controlling the COVID-19 outbreak. The purpose of this study was to set up a Sanger-based platform for massive SARS-CoV-2 variant tracking in laboratories in low-resource settings. METHODS: We used nested RT-PCR assay, Sanger sequencing and lineage assignment for 930-bp of the SARS-CoV-2 spike gene, which harbors specific variants of concern (VOCs) mutations. We set up our platform by comparing its results with whole genome sequencing (WGS) data on 137 SARS-CoV-2 positive samples. Then, we applied it on 1028 samples from March-September 2021. RESULTS: In total, 125 out of 137 samples showed 91.24% concordance in mutation detection. In lineage assignment, 123 out of 137 samples demonstrated 89.78% concordance, 65 of which were assigned as VOCs and showed 100% concordance. Of 1028 samples screened by our in-house method, 78 distinct mutations were detected. The most common mutations were: S:D614G (21.91%), S:P681R (12.19%), S:L452R (12.15%), S:T478K (12.15%), S:N501Y (8.91%), S:A570D (8.89%), S:P681H (8.89%), S:T716I (8.74%), S:L699I (3.50%) and S:S477N (0.28%). Of 1028 samples, 980 were attributed as VOCs, which include the Delta (B.1.617.2) and Alpha (B.1.1.7) variants. CONCLUSION: Our proposed in-house Sanger-based assay for SARS-CoV-2 lineage assignment is an accessible strategy in countries with poor infrastructure facilities. It can be applied in the rapid tracking of SARS-CoV-2 VOCs in the SARS-CoV-2 pandemic.

Bi-allelic loss of function variant in the NRCAM gene is associated with motor-predominant axonal polyneuropathy; the second report.

BACKGROUND: The role of biallelic variants in the NRCAM gene underlying a neurodevelopmental disorder has been defined recently. The phenotype is mainly recognized by varying severity of global developmental delay/intellectual disability, hypotonia, spasticity, and peripheral neuropathy. METHODS: Here, we describe a patient with an initial diagnosis of motor-predominant axonal polyneuropathy or a form of distal SMA. Whole-exome sequencing (WES), in parallel with WES-based CNV detection and assessment of homozygosity runs, was performed to identify this patient's possible genetic cause. RESULTS: Whole exome sequencing revealed a homozygous variant, c.73C > T (p.Gln25*), in the NRCAM gene, while the patient manifests a mild range of phenotypes compared to NRCAM-related disorder. He presented only motor-predominant axonal polyneuropathy with no other signs of central nervous system involvement. CONCLUSIONS: This study is the second report of an association between biallelic NRCAM gene variants and a Mendelian disorder. The obtained clinical data, together with the molecular findings in this patient, expands the clinical and molecular spectrum of NRCAM-related disorder and highlights its phenotypic complexity. Although patients with loss of function variants in this gene have previously presented severe clinical features, we show that type of the pathogenic variant does not necessarily determine the severity of this phenotype.

Emerging Epidemiological Data on Rare Intellectual Disability Syndromes from Analyzing the Data of a Large Iranian Cohort.

BACKGROUND: Intellectual disability (ID) is a genetically heterogeneous condition, and so far, 1679 human genes have been identified for this phenotype. Countries with a high rate of parental consanguinity, such as Iran, provide an excellent opportunity to identify the remaining novel ID genes, especially those with an autosomal recessive (AR) mode of inheritance. This study aimed to investigate the most prevalent ID genes identified via next-generation sequencing (NGS) in a large ID cohort at the Genetics Research Center (GRC) of the University of Social Welfare and Rehabilitation Sciences. METHODS: First, we surveyed the epidemiological data of 619 of 1295 families in our ID cohort, who referred to the Genetics Research Center from all over the country between 2004 and 2021 for genetic investigation via the NGS pipeline. We then compared our data with those of several prominent studies conducted in consanguineous countries. Data analysis, including cohort data extraction, categorization, and comparison, was performed using the R program version 4.1.2. RESULTS: We categorized the most common ID genes that were mutated in more than two families into 17 categories. The most common syndromic ID in our cohort was AP4 deficiency syndrome, and the most common non-syndromic autosomal recessive intellectual disability (ARID) gene was ASPM. We identified two unrelated families for the 36 ID genes. We found 14 genes in common between our cohort and the Arab and Pakistani groups, of which three genes (AP4M1, AP4S1, and ADGRG1) were repeated more than once. CONCLUSION: To date, there has been no comprehensive targeted NGS platform for the detection of ID genes in our country. Due to the large sample size of our study, our data may provide the initial step toward designing an indigenously targeted NGS platform for the diagnosis of ID, especially common ARID in our population.

Expanding the Molecular Spectrum of HK1-Related Charcot-Marie-Tooth Disease, Type 4G; the First Report in Iran.

Charcot-Marie-Tooth disease type 4G (CMT4G) was first reported in Balkan Gypsies as a myelinopathy starting with progressive distal lower limb weakness, followed by upper limb involvement and prominent distal sensory impairment later in the patient's life. So far, CMT4G has been only reported in European Roma communities with two founder homozygous variants; g.9712G>C and g.11027G>A, located in the 5'-UTR of the HK1 gene. Here, we present the first Iranian CMT4G patient manifesting progressive distal lower limb weakness from 11 years of age and diagnosed with chronic demyelinating sensorimotor polyneuropathy. Whole-exome sequencing for this patient revealed a homozygous c.19C>T (p. Arg7*) variant in the HK1 gene. This report expands the mutational spectrum of the HK1-related CMT disorder and provides supporting evidence for the observation of CMT4G outside the Roma population. Interestingly, the same Arg7* variant is recently observed in another unrelated Pakistani CMT patient, proposing a possible prevalence of this variant in the Middle Eastern populations.

ZBTB11 dysfunction: spectrum of brain abnormalities, biochemical signature and cellular consequences.

Bi-allelic pathogenic variants in ZBTB11 have been associated with intellectual developmental disorder, autosomal recessive 69 (MRT69; OMIM 618383). We report five patients from three families with novel, bi-allelic variants in ZBTB11. We have expanded the clinical phenotype of MRT69, documenting varied severity of atrophy affecting different brain regions and described combined malonic and methylmalonic aciduria as a biochemical manifestation. As ZBTB11 encodes for a transcriptional regulator, we performeded chromatin immunoprecipitation-sequencing targeting ZBTB11 in fibroblasts from patients and controls. Chromatin immunoprecipitation-sequencing revealed binding of wild-type ZBTB11 to promoters in 238 genes, among which genes encoding proteins involved in mitochondrial functions and RNA processing are over-represented. Mutated ZBTB11 showed reduced binding to 61 of the targeted genes, indicating that the variants act as loss of function. Most of these genes are related to mitochondrial functions. Transcriptome analysis of the patient fibroblasts revealed dysregulation of mitochondrial functions. In addition, we uncovered that reduced binding of the mutated ZBTB11 to ACSF3 leads to decreased ACSF3 transcript level, explaining combined malonic and methylmalonic aciduria. Collectively, these results expand the clinical spectrum of ZBTB11-related neurological disease and give insight into the pathophysiology in which the dysfunctional ZBTB11 affect mitochondrial functions and RNA processing contributing to the neurological and biochemical phenotypes.

Disease Waves of SARS-CoV-2 in Iran Closely Mirror Global Pandemic Trends.

BACKGROUND: Complete SARS-CoV-2 genome sequencing in the early phase of the outbreak in Iran showed two independent viral entries. Subsequently, as part of a genome surveillance project, we aimed to characterize the genetic diversity of SARS-CoV-2 in Iran over one year after emerging. METHODS: We provided 319 SARS-CoV-2 whole-genome sequences used to monitor circulating lineages in March 2020-May 2021 time interval. RESULTS: The temporal dynamics of major SARS-CoV-2 clades/lineages circulating in Iran is comparable to the global perspective and represent the 19A clade (B.4) dominating the first disease wave, followed by 20A (B.1.36), 20B (B.1.1.413), 20I (B.1.1.7), leading the second, third and fourth waves, respectively. We observed a mixture of circulating B.1.36, B.1.1.413, B.1.1.7 lineages in winter 2021, paralleled in a fading manner for B.1.36/B.1.1.413 and a growing rise for B.1.1.7, prompting the fourth outbreak. Entry of the Delta variant, leading to the fifth disease wave in summer 2021, was detected in April 2021. This study highlights three lineages as hallmarks of the SARS-CoV-2 outbreak in Iran; B4, dominating early periods of the epidemic, B.1.1.413 (B.1.1 with the combination of [D138Y-S477N-D614G] spike mutations) as a characterizing lineage in Iran, and the co-occurrence of [I100T-L699I] spike mutations in half of B.1.1.7 sequences mediating the fourth peak. It also designates the renowned combination of G and GR clades' mutations as the top recurrent mutations. CONCLUSION: In brief, we provided a real-time and comprehensive picture of the SARS-CoV-2 genetic diversity in Iran and shed light on the SARS-CoV-2 transmission and circulation on the regional scale.

SARS-CoV-2 outbreak in Iran: The dynamics of the epidemic and evidence on two independent introductions.

The SARS-CoV-2 virus has been rapidly spreading globally since December 2019, triggering a pandemic, soon after its emergence. While Iran was among the first countries confronted with rapid spread of virus in February 2020, no real-time SARS-CoV-2 whole-genome tracking in early phase of outbreak was performed in the country. To address this issue, we provided 50 whole-genome sequences of viral isolates ascertained from different geographical locations in Iran during March-July 2020. The corresponding analysis on origins, transmission dynamics and genetic diversity of SARS-CoV-2 virus, represented at least two introductions of the virus into the country, constructing two major clusters defined as B.4 and B.1*. The first entry of the virus might have occurred around very late 2019/early 2020, as suggested by the time to the most recent common ancestor, followed by a rapid community transmission that led to dominancy of B.4 lineage in early epidemic till the end of June. Gradually, reduction in dominancy of B.4 occurred possibly as a result of other entries of the virus, followed by surge of B.1* lineages, as of mid-May. Remarkably, variation tracking of the virus indicated the increase in frequency of D614G mutation, along with B.1* lineages, which showed continuity till October 2020. The increase in frequency of D614G mutation and B.1* lineages from mid-May onwards predicts a rapid viral transmission that may push the country into a critical health situation followed by a considerable change in composition of viral lineages circulating in the country.

CEP104 and CEP290; Genes with Ciliary Functions Cause Intellectual Disability in Multiple Families.

BACKGROUND: Neurodevelopmental and intellectual impairments are extremely heterogeneous disorders caused by a diverse variety of genes involved in different molecular pathways and networks. Genetic alterations in cilia, highly-conserved organelles with sensorineural and signal transduction roles can compromise their proper functions and lead to so-called "ciliopathies" featuring intellectual disability (ID) or neurodevelopmental disorders as frequent clinical manifestations. Here, we report several Iranian families affected with ID and other ciliopathy-associated features carrying known and novel variants in two ciliary genes; CEP104 and CEP290. METHODS: Whole exome and Targeted exome sequencing were carried out on affected individuals. Lymphoblastoid cell lines (LCLs) derived from the members of affected families were established for two families carrying CEP104 mutations. RNA and protein expression studies were carried out on these cells using qPCR and Western blot, respectively. RESULTS: A novel homozygous variant; NM_025114.3:c.7341_7344dupACTT p.(Ser2449Thrfs*8) and four previously reported homozygous variants; NM_025114.3:c.322C>T p.(Arg108*), NM_025114.3:c.4393C>T p.(Arg1465*), NM_025114.3:c.5668G>T p.(Gly1890*) and NM_025114.3:c.1666dupA p.(Ile556Asnfs*20) were identified in CEP290. In two other families, two novel homozygous variants; NM_014704:c.2356_2357insTT p.(Cys786Phefs*11) and NM_014704:c.1901_1902insT p.(Leu634Phefs*33) were identified in CEP104, another ciliary gene. qPCR and Western blot analyses showed significantly lower levels of CEP104 transcripts and protein in patients compared to heterozygous or normal family members. CONCLUSION: We emphasize on the clinical variability and pleiotropic phenotypes due to variants of these genes. In conclusion, our findings support the pivotal role of these genes resulting in cognitive and neurodevelopmental features.

Biallelic mutations in the death domain of PIDD1 impair caspase-2 activation and are associated with intellectual disability.

PIDD1 encodes p53-Induced Death Domain protein 1, which acts as a sensor surveilling centrosome numbers and p53 activity in mammalian cells. Early results also suggest a role in DNA damage response where PIDD1 may act as a cell-fate switch, through interaction with RIP1 and NEMO/IKKg, activating NF-κB signaling for survival, or as an apoptosis-inducing protein by activating caspase-2. Biallelic truncating mutations in CRADD-the protein bridging PIDD1 and caspase-2-have been reported in intellectual disability (ID), and in a form of lissencephaly. Here, we identified five families with ID from Iran, Pakistan, and India, with four different biallelic mutations in PIDD1, all disrupting the Death Domain (DD), through which PIDD1 interacts with CRADD or RIP1. Nonsense mutations Gln863* and Arg637* directly disrupt the DD, as does a missense mutation, Arg815Trp. A homozygous splice mutation in the fifth family is predicted to disrupt splicing upstream of the DD, as confirmed using an exon trap. In HEK293 cells, we show that both Gln863* and Arg815Trp mutants fail to co-localize with CRADD, leading to its aggregation and mis-localization, and fail to co-precipitate CRADD. Using genome-edited cell lines, we show that these three PIDD1 mutations all cause loss of PIDDosome function. Pidd1 null mice show decreased anxiety, but no motor abnormalities. Together this indicates that PIDD1 mutations in humans may cause ID (and possibly lissencephaly) either through gain of function or secondarily, due to altered scaffolding properties, while complete loss of PIDD1, as modeled in mice, may be well tolerated or is compensated for.

Comprehensive genotype-phenotype correlation in AP-4 deficiency syndrome; Adding data from a large cohort of Iranian patients.

Mutations in adaptor protein complex-4 (AP-4) genes have first been identified in 2009, causing a phenotype termed as AP-4 deficiency syndrome. Since then several patients with overlapping phenotypes, comprised of intellectual disability (ID) and spastic tetraplegia have been reported. To delineate the genotype-phenotype correlation of the AP-4 deficiency syndrome, we add the data from 30 affected individuals from 12 out of 640 Iranian families with ID in whom we detected disease-causing variants in AP-4 complex subunits, using next-generation sequencing. Furthermore, by comparing genotype-phenotype findings of those affected individuals with previously reported patients, we further refine the genotype-phenotype correlation in this syndrome. The most frequent reported clinical findings in the 101 cases consist of ID and/or global developmental delay (97%), speech disorders (92.1%), inability to walk (90.1%), spasticity (77.2%), and microcephaly (75.2%). Spastic tetraplegia has been reported in 72.3% of the investigated patients. The major brain imaging findings are abnormal corpus callosum morphology (63.4%) followed by ventriculomegaly (44.5%). Our result might suggest the AP-4 deficiency syndrome as a major differential diagnostic for unknown hereditary neurodegenerative disorders.

Molecular Diagnosis of Hereditary Neuropathies by Whole Exome Sequencing and Expanding the Phenotype Spectrum.

BACKGROUND: Inherited peripheral neuropathies (IPNs) are a group of neuropathies affecting peripheral motor and sensory neurons. Charcot-Marie-Tooth (CMT) disease is the most common disease in this group. With recent advances in next-generation sequencing (NGS) technologies, more than 100 genes have been implicated for different types of CMT and other clinically and genetically inherited neuropathies. There are also a number of genes where neuropathy is a major feature of the disease such as spinocerebellar ataxia (SCA) and hereditary spastic paraplegia (HSP). We aimed to determine the genetic causes underlying IPNs in Iranian families. METHODS: We performed whole exome sequencing (WES) for 58 PMP22 deletion-/duplication-negative unrelated Iranian patients with a spectrum of phenotypes and with a preliminary diagnosis of hereditary neuropathies. RESULTS: Twenty-seven (46.6%) of the cases were genetically diagnosed with pathogenic or likely pathogenic variants. In this study, we identified genetically strong variants within genes not previously linked to any established disease phenotype in five (8.6%) patients. CONCLUSION: Our results highlight the advantage of using WES for genetic diagnosis in highly heterogeneous diseases such as IPNs. Moreover, functional analysis is required for novel and uncertain variants.

Whole genome sequencing identifies a duplicated region encompassing Xq13.2q13.3 in a large Iranian family with intellectual disability.

BACKGROUND: The X chromosome has historically been one of the most thoroughly investigated chromosomes regarding intellectual disability (ID), whose etiology is attributed to many factors including copy number variations (CNVs). Duplications of the long arm of the X chromosome have been reported in patients with ID, short stature, facial anomalies, and in many cases hypoplastic genitalia and/or behavioral abnormalities. METHODS: Here, we report on a large Iranian family with X-linked ID caused by a duplication on the X chromosome identified by whole genome sequencing in combination with linkage analysis. RESULTS: Seven affected males in different branches of the family presented with ID, short stature, seizures, facial anomalies, behavioral abnormalities (aggressiveness, self-injury, anxiety, impaired social interactions, and shyness), speech impairment, and micropenis. The duplication of the region Xq13.2q13.3, which is ~1.8 Mb in size, includes seven protein-coding OMIM genes. Three of these genes, namely SLC16A2, RLIM, and NEXMIF, if impaired, can lead to syndromes presenting with ID. Of note, this duplicated region was located within a linkage interval with a LOD score >3. CONCLUSION: Our report indicates that CNVs should be considered in multi-affected families where no candidate gene defect has been identified in sequencing data analysis.

Brief Report of Variants Detected in Hereditary Hearing Loss Cases in Iran over a 3-Year Period.

BACKGROUND: Diagnosis of hereditary hearing loss (HHL) as a heterogeneous disorder is very important especially in countries with high rates of consanguinity where the autosomal recessive pattern of inheritance is prevalent. Techniques such as next-generation sequencing, a comprehensive genetic test using targeted genomic enrichment and massively parallel sequencing (TGE + MPS), have made the diagnosis more cost-effective. The aim of this study was to determine HHL variants with comprehensive genetic testing in our country. METHODS: Fifty GJB2 negative individuals with HHL were referred to the Kariminejad-Najmabadi Pathology and Genetics Center, Tehran, one of the reference diagnostic genetic laboratories in Iran, during a 3-year period between 2014 and 2017. They were screened with the OtoSCOPE test, the targeted genomic enrichment and massively parallel sequencing (TGE + MPS) platform after a detailed history had been taken along with clinical evaluation. RESULTS: Among 32 out of 50 GJB2 negative patients (64%), 34 known pathogenic and novel variants were detected of which 16 (47%) were novel, identified in 10 genes of which the most prevalent were CDH23, MYO7A and MYO15A. CONCLUSION: These results provide a foundation from which to make appropriate recommendations for the use of comprehensive genetic testing in the evaluation of Iranian patients with hereditary hearing loss.

Distinct genetic variation and heterogeneity of the Iranian population.

Iran, despite its size, geographic location and past cultural influence, has largely been a blind spot for human population genetic studies. With only sparse genetic information on the Iranian population available, we pursued its genome-wide and geographic characterization based on 1021 samples from eleven ethnic groups. We show that Iranians, while close to neighboring populations, present distinct genetic variation consistent with long-standing genetic continuity, harbor high heterogeneity and different levels of consanguinity, fall apart into a cluster of similar groups and several admixed ones and have experienced numerous language adoption events in the past. Our findings render Iran an important source for human genetic variation in Western and Central Asia, will guide adequate study sampling and assist the interpretation of putative disease-implicated genetic variation. Given Iran's internal genetic heterogeneity, future studies will have to consider ethnic affiliations and possible admixture.

Iranome: A catalog of genomic variations in the Iranian population.

Considering the application of human genome variation databases in precision medicine, population-specific genome projects are continuously being developed. However, the Middle Eastern population is underrepresented in current databases. Accordingly, we established Iranome database (www.iranome.com) by performing whole exome sequencing on 800 individuals from eight major Iranian ethnic groups representing the second largest population of Middle East. We identified 1,575,702 variants of which 308,311 were novel (19.6%). Also, by presenting higher frequency for 37,384 novel or known rare variants, Iranome database can improve the power of molecular diagnosis. Moreover, attainable clinical information makes this database a good resource for classifying pathogenicity of rare variants. Principal components analysis indicated that, apart from Iranian-Baluchs, Iranian-Turkmen, and Iranian-Persian Gulf Islanders, who form their own clusters, rest of the population were genetically linked, forming a super-population. Furthermore, only 0.6% of novel variants showed counterparts in "Greater Middle East Variome Project", emphasizing the value of Iranome at national level by releasing a comprehensive catalog of Iranian genomic variations and also filling another gap in the catalog of human genome variations at international level. We introduce Iranome as a resource which may also be applicable in other countries located in neighboring regions historically called Greater Iran (Persia).